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Abstract: The recently proposed hard-pion chiral perturbation theory predicts that the 
leading chiral logarithms factorize with respect to the energy dependence in the chiral limit. 
This claim has been successfully tested in the pion form factors up to two loops in chiral 
perturbation theory. In the present paper we explain this factorization property at two loops 
and even show that it is valid to all orders for a subclass of diagrams. We also demonstrate 
that factorization is violated starting at three loops. 
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1 Introduction 

Chiral perturbation theory (%PT) [1, 2] allows one to calculate the momentum and light quark 
mass dependence of hadronic processes at low energy. It can also be applied to processes in 
which hard scales are involved, e.g. when baryon or heavy quarks interact with pions. In 
these cases, however, the range of applicability of the theory is restricted to the kinematical 
regions in which the momenta of the pions are soft. In ref. [3] Flynn and Sachrajda proposed 
an interesting argument which, after reinterpreting the effective Lagrangian describing the 
decay of a heavy meson (the K in that case) into a pion and a lepton- neutrino pair, leads to a 
prediction of the coefficient of the leading chiral logarithm in decays even for kinematical 
configurations where the pion is not soft. Subsequently, in a series of papers by Bijnens 
and collaborators [4-7] it has been claimed that the calculation of such a chiral logarithm is 
possible in more general processes in which the pion is hard. This approach has been referred 
to as hard pion chiral perturbation theory (H7rxPT). 

A particularly clear example of a prediction in HvrxPT has been provided in ref. [6] 
where it has been shown that the scalar and vector form factors of the pion, which have been 
calculated to two loops in SU(2) %PT [8], in the limit M%/s -C 1 factorize: 

F ViS (s, Ml) = F v>s (s) (1 + <*vA*) L ) + 0{ M 2 ) ■ ( L1 ) 
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Here s is the momentum transfer squared and L stands for the leading chiral logarithm, 
defined as 1 

M 2 . M 2 . . 

M 2 is proportional to the average up and down quark masses m (M 2 = 2Bm) and F is the 
pion decay constant in the chiral limit: 

M 2 = M 2 + C(M 4 ) , F W = F + 0{M 2 ) . (1.3) 

Fv,s( s ) denote the form factors for vanishing u and d quark masses 2 . Bijnens and Jemos 
provide arguments in support of the validity of this factorization property to all orders [6]. 
On the other hand a precise formulation of this factorization property is still lacking and it 
is not clear under what circumstances and for what quantities it holds. 

The aim of this paper is to start providing an answer to these questions. In particular 
we asked ourselves whether factorization of chiral logs for asymptotically small values of the 
quark masses could be exact, or in other words, a property of QCD. If the latter is true, 
then this property must emerge also in any effective theory of QCD provided this can deal 
with the chiral limit. This is the case of m the limit M 2 <C s <C A 2 , and indeed the 
observed factorization in the form factors at two-loop level points in this direction. In the 
present paper we will analyze the pion form factors in beyond two loops. We find that 
the most effective way to attack this problem is to study the dispersion relation for the form 
factors and apply the chiral counting to this. As we will show, this allows us to formulate a 
recursive argument and to establish factorization for a subclass of diagrams to all orders in 
the chiral counting. 

We also find, however, that inelastic contributions (three-loop diagrams with a four-pion 
intermediate state) also generate chiral logs of the order of interest to us and that these violate 
the factorization property. These diagrams do not fall in the subclass of diagrams which are 
treated in our recursive argument mentioned above. We conclude that factorization of chiral 
logs cannot be an exact property of QCD for the pion form factors. If we believe that 
HzrxPT is valid even for (asymptotically) large values of s, then this conclusion would be 
implied by the available results on the meson form factors in the limit of asymptotically large 
energy. In this limit the pion form factors can be analyzed in the context of perturbative 
QCD factorization, as shown by Brodsky-Lepage [9]. In their formula the quark masses are 
neglected, but as discussed by Chen and Stewart [10], it is possible to study the deviation 
from the chiral limit and in particular to determine the coefficient of the leading chiral log. 
This chiral log factorizes with respect to the leading Brodsky-Lepage term (but there is no 
reason to assume that it could factorize for the subleading terms). However, the coefficient 
of the leading chiral log is different from what is obtained in the low-energy limit. 

1 Writing lnAf 2 /s as lnM 2 //i 2 + ln/i 2 /s, one can then equivalently define L as /^-dependent as done in 
refs. [3, 6] and the second term goes into the 0(M 2 ) contribution in eq. (1.1). 

2 In this paper quantities in the SU(2) chiral limit are denoted by a bar, X — liniAf^-xD X. 
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The paper is organized as follows: in sec. 2 we introduce our notation and briefly discuss 
the dispersive representation of the pion form factors. In sec. 3 we show how, after applying 
the chiral counting, one can determine the leading chiral log on the basis of the dispersive 
representation. Considering only elastic contributions (two-pion intermediate states) we are 
able to provide a recursive argument and prove factorization of the leading chiral log to all 
orders. In sec. 4 we discuss the leading inelastic contributions, calculate the chiral log coming 
from these and show that these violate factorization. In sec. 5 we briefly discuss the papers 
of Brodsky-Lepage [9] and Chen-Stewart [10] and show that these also lead to the same 
conclusion, namely that factorization of the leading chiral log cannot be exact. 



2 Dispersive representation of the pion form factors 

We consider the vector and scalar pion form factors, respectively Fy(s) and Fg(s) and denote 
them by the common symbol F(s) (unless it is necessary to distinguish them). We normalize 
the scalar form factor such that 3 Fs(0) = 1 — for the vector form factor this condition follows 
from current conservation. Note that this normalization condition (whether automatic or 
imposed) is inessential to our argument. Both these form factors are analytic functions in the 
cut plane [4M^,oo) and satisfy a once subtracted dispersion relation of the form 



s f°° , , ImF(s') , , 

F(s) = l + - ds' / \ . 2.1 
7T Jam* s (« - s) 



The subtraction constant is the value of the form factor at zero momentum transfer which is 
equal to one as mentioned above. In the elastic region unitarity relates the imaginary part 
of the form factor to the form factor itself and the hit partial wave with the same quantum 
numbers: 

/ 4M 2 

lm F(s) = a(s)F(s)t*(s) , ff ( s ) = Wl--^. (2.2) 

When s gets larger than the inelastic threshold, additional cuts involving several intermediate 
pions contribute to the discontinuity. In the subsequent discussion we will consider the form 
factor which is obtained by solving the dispersion relation with imaginary parts arising only 
from two-pion intermediate states. We will be able to make a statement about this con- 
tribution, provided we consider only two-pion intermediate states also for the tttt scattering 
amplitude entering the unitarity relation (2.2). We call this part the elastic contribution, and 
the rest the "inelastic" part. 

F(s) = F el (s) + F inel (s) . (2.3) 

We stress that this definition is a diagrammatic one. Its precise meaning should become more 
clear in the following section, in which we will consider only the elastic contribution. The 
first diagrams contributing to i*i ne i(s), which start at order p 6 , are discussed in sec. 4. 



3 Notice that, consistent with this normalization, we define leading order as O(p ), next-to-leading order as 
0(p 2 ) and so on. The reader should be aware of the mismatch with the usual counting for the Lagrangian and 
the tvty scattering amplitude discussed below, where leading order is 0(p 2 ). 
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In the following we will apply the chiral counting to the dispersion relation [11, 12]. 
It is well-known 4 that if one does so one gets a one-to-one correspondence between direct 
xPT calculations and the various chiral orders in the dispersion relation: in other words 
xPT provides a perturbative solution to the dispersion relation. This will be the basis of 
our analysis 5 . Indeed, since t(s) starts at 0(p 2 ), knowing the form factor and the tvk partial 
wave up to a given chiral order means to know Im F(s) at one order higher. If one is then 
able to perform the dispersive integral, the corresponding real part of the form factor is also 
obtained. This is a powerful tool which allows us to argue recursively. 

The chiral representation contains of course more information than what one gets from 
the dispersion relation as it also gives the quark mass dependence of subtraction constants. 
Indeed whenever necessary we will use xPT to determine chiral logs in subtraction constants. 
This will be needed in the ttit scattering amplitude, but not for the form factors since with 
the normalization condition we chose for the form factors (2.1), the subtraction constant is 
equal to one and cannot contain any logs. 

3 The elastic contribution to F(s) 

We will now address the question, how leading chiral logarithms can arise from the elastic 
contribution to the dispersive integral. There are two possible mechanisms to generate terms 
containing L: either one starts from an integrand which does not contain a log of the pion 
mass and this is produced by the integration over s', or the integrand itself contains a chiral 
log and the dispersive integral determines what function of s multiplies it. We will now 
discuss these two mechanisms in turn. We stress that, as mentioned in the introduction, 
we are considering the form factors in the limit M 2 <C s and are only interested in terms 
proportional to L: terms of 0(M 2 ) without logarithmic enhancement or terms of 0(M 2 L) 
are all beyond the accuracy we aim at. 

To avoid clutter in the notation in this section we will drop the subscript "el", F(s) = 
F e i(s), as this is the only contribution discussed here. 

3.1 Leading chiral logarithms from the dispersive integration 

The first possible mechanism is that the leading chiral logarithms are generated by the dis- 
persive integral, i.e. they are produced at the lower integration boundary s' ~ 4M^, which 
goes to zero in the chiral limit. In order to investigate this, we must analyze the behavior of 
the integrand in the regime s 1 ~ M 2 <C s: we can therefore make use of the standard chiral 
expansion for the discontinuity of the form factor (2.2). Expanding both the form factor and 
the 7T7T amplitude we can write 

F(s,M 2 ) = l + 0(s) + 0{M 2 ) and t(s, M 2 ) = c x M 2 + c 2 s + 0(s 2 ) + 0{M 4 ) . (3.1) 
4 For a pedagogical discussion see [13]. 

5 Early calculations of the chiral logs in chiral perturbation theory have been made following a similar 
approach [14, 15]. 
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Plugging eq. (3.1) into eq. (2.1) we obtain 

F(s) = 1 + - ds' ° { f ] U M 2 + c 2 s' + 0(p 4 )) . (3.2) 

The three terms in brackets in the integrand in eq. (3.2) generate three types of integrals: 
the first (the one proportional to c\) is the well-known loop function J{s) which has the 
following expansion in M 2 /s < 1, 

at 00 ,, o(s') — . , 1 |" M 2 2M 2 ( , M 2 \ m f M 4 \] . . 

~ ds ,t ) \ = 16ttJ(s = - 1 + ln + 1-ln + O I ) . 3.3 

7T Jam* s'(s'-s) tt [ -s s \ -s J V s /. 

The remaining two integrals are UV divergent. However, since we are interested in their 
behaviour close to the lower integration boundary, we can introduce an M^-independent cut- 
off A 2 , which allows us to interchange integration and expansion for small M. The second 
type of integral (the one proportional to C2) is then given by 



[ A ds'j^r = d 1 ( s ,A 2 )--M 2 ln— + 0(M 2 ) (3.4) 



and the third by 



2 



, ds' ,. ' . x 0(p 4 ) = d 2 (s, A 2 ) + 0(M 2 ) . (3.5) 
n Jam* s'(s'-s) 

With 0(p ) we indicate terms proportional to s' 2 , s'M 2 and M , which are all small in the 
region close to the lower integration boundary. It is easy to see that this further suppression 
does not allow terms proportional to L (without further powers of M 2 ) to be generated. We 
conclude that the dispersive integral can generate leading chiral logs only from the leading 
chiral contribution to the integrand. The functions di(s,A 2 ) and A 2 ) are quark mass 
independent and contribute to F(s) together with the low-energy constants (LECs) which 
cancel the A-dependence. 

Putting all pieces together we conclude that the chiral log generated by the dispersive 
integration is given by 

16vrF 2 (ci - 2c 2 )L . (3.6) 

The constants c\ and c 2 are related to the leading chiral contributions to the tttx scattering 
lengths and effective ranges characterizing the threshold expansion 

Ret'is) = q 2£ [4 + b\q 2 + 0(g 4 )) , (3.7) 

where q 2 = s/A — M 2 and I and I denote isospin and angular momentum, respectively. For the 
scalar form factor the relevant parameters are ci = a[j/M 2 — b^+0{M 2 ) and c 2 = 6g/4+0(M 2 ) 
where [2] 

7M 2 1 
°° = 32^ + ° (M4) ' b °° = 4^ + 0{M2) ' (3 ' 8) 
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leading to 



a s = 16vrF 2 (ci - 2 c 2 ) = 16ttF 2 lim ( - 



3bg 



mCo 2 ' - ' (3 ' 9) 

which reproduces the known result [6, 8]. For the vector form factor we must instead use 
the parameters: c\ = — (a\)\ 3 = — l/(247rF 2 ), and c 2 = (°i/4)| 2 = l/(967rF 2 ), which 
leads to 



a v = 16ttF 2 (ci - 2c 2 ) 



16nF 2 lim 

M 2 -s>0 



3a} 



-1 



(3.10) 



also in agreement with the explicit calculation in refs. [6, 8]. 

For the subsequent discussion it is useful to determine the leading order of F(s) in the 
chiral expansion in powers of s (see also [16]). In order to do this, we have to find a function 
which is analytic in the cut plane [0, oo) and with the following imaginary part along the cut 



lim 

M->0 



)F^(s)t^( 



c 2 s 



Such a function is easily found: 



F (s) = s — In — 

7T — S 



(3.11) 



(3.12) 



with A 2 an unknown scale. The explicit expressions in the case of the scalar and vector form 
factors read 



16tt 2 F 2 
s 

16vr 2 F 2 



1 + In — + i-w + 16ttX(u) 
s 



1, 



?.7T 



+ - In — + — - 16k 2 eJu) 
18 6 s 6 6 



(3.13) 



in agreement with ref. [6]. The coefficients of the logarithms are indeed correctly repro- 
duced by substituting c 2 = (^o/4)| A/2 ^ = l/(167rF 2 ) for the scalar form factor and c 2 = 
(aJ/4)[ 2 = 1/(96ttF 2 ) for the vector one. 

3.2 Leading chiral logarithms from the integrand 

We shall now discuss whether and how leading chiral logarithms can be generated at higher 
chiral orders — more specifically, we are interested in terms proportional to s n ~ 1 L at order 
p 2n . The discussion above made it clear that the behavior of the integrand around the 
lower limit of integration may only generate a chiral logarithm at 0(p 2 ). There is a second 
mechanism, however, by which chiral logarithms may arise from the dispersive integrals at 
higher orders, namely if the integrand itself contains a chiral logarithm. 

Let us consider the dispersion relation at 0(p A ) (i.e. the contributions to the form factor 
at the two- loop level). At this order the integrand has this form 



Imf( 4 )(s) = a{s) \t^*(s) + F ( - 2 \s)t^(s) 



(3.14) 
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and each of the terms may contain chiral logarithms. We consider first the latter one: 

is the tree-level contribution to the tttt scattering amplitude and contains no chiral logarithms, 

whereas F^ 2 \s) does, as we have seen above. Expanding in M 2 /s, we can write it as 

F^ 2 \s) tW{a) = (f (2 \s) + aLj t {2 \s) + 0(M 2 ) . (3.15) 

We can similarly expand t^\s), 

t { - A Xs)=t {A \s) + PsL + 0{M 2 ) . (3.16) 

It is easy to realize that the term proportional to j3 would destroy factorization: this vanishes, 
however, as shown in app. A. We therefore conclude that the only term containing L in 
eq. (3.14) has as coefficient exactly the absorptive part of the form factor at one chiral order 
lower times a. As we have discussed above, the solution of the corresponding dispersion 
relation reads 2 

as Sl\ n ±x (3.17) 

7T —S 

with Ax an unknown energy scale. We argue, however, that this has to coincide with A2 
introduced in eq. (3.12). From the point of view of xPT both of these scales are related to 
LECs in the 0(p 4 ) La grangian: A2 in particular is the one which appears in the form factor 
at this order. Since we have just showed that one cannot generate a chiral logarithm by 
integrating over a local contribution of 0(p 4 ) to the tttt scattering amplitude (which is the 
only other vertex in the dispersion relation), we have to conclude that Ax = A2. We conclude 
that at two loops we can write the contribution to the form factor as 

F(s) = (l + F i2) (s)) (l + aL) + F W (s) + 0(M 2 ) + 0(p 6 ) . (3.18) 



i.e. in a factorized form, as predicted by H7rxPT. 

For the elastic contribution to F(s), the same reasoning can be repeated in exactly the 
same way order by order. At every new step the terms threatening factorization are the 
contributions to ImF' 2 "' arising from the tttt scattering amplitude at the same order. A 
chiral logarithm of the form s n ~ l L in t^ 2 ^ would destroy factorization. In app. A we show 
that these terms are absent. However, at order p 6 , the four-pion intermediate states contribute 
to Im F(s) and we will now show that these yield leading chiral logarithms and are responsible 
for the breakdown of factorization, leading to the three-loop result 

F(s) = (l + F {2) (s)+F { *\s)) (l + aL) + a iQcl (s)L + F i6 \s) + 0(M 2 ) + 0(p 8 ) . (3.19) 



Before closing the section we comment on the form factor in the chiral limit: at 0(p 4 ) 
this is given by the solution of the dispersion relation with discontinuity 



Im F (i \s)= t w *( S ) + F w ( S )F ; (s) . (3.20) 
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The form factor F^\s) can be derived from ref. [8] and given in explicit analytic form. The 
expression for the scalar form factor, which will be used for the numerical analysis in sec. 4, 
e.g. reads 

"" ""' ' (3-21) 



where N := 16vr 2 , £ r a := 1/3(11*1 + 7^) +£ r 4 and Z := 1/JV (hi(s//i 2 ) - ivr - l). The r£ 3 LEC 
stems from the 0(p 6 ) chiral Lagrangian. 

4 The contribution from inelastic channels 

In the previous section we have shown that the leading chiral logarithms originating from 
the elastic part factorize. We are now going to show that the 0(p G ) inelastic contribution 
leads also to terms proportional to s 2 L, i.e. ai ne i(s) ^ in eq. (3.19), which implies that 
the factorization hypothesis at the basis of HvrxPT is not valid to all orders. We shall now 
present the details of our calculation. 

The lowest order inelastic contribution to F(s) is given by three-loop diagrams with four 
intermediate pions. We evaluated them by means of the following dispersion relation with 
the lower integration boundary given by the four-pion threshold: 

8 f°° , ImF incl ( S ') 
^inel(s) = - / ds — — . (4.1) 

Chiral logarithms are either produced by the dispersive integration or are contained in the 
integrand. In the chiral counting Im.Fj ne i is of C(p 6 ) and contains a four-particle phase space 
factor which gives a strong suppression near threshold ( <&i{s) ~ (s — 16M 2 ) 7 / 2 ). Arguing 
similarly to how we did in sec. 3.1 we conclude that leading chiral logs cannot be generated 
by the dispersive integration. Therefore we need to concentrate on the integrand. 
Due to unitarity, the imaginary part of the form factor has this form 

ImF incl (s) = i J d<S> 4 (s;pi,p2,P3,P4) F 4ir ■ T* n , (4.2) 

where d<3?4 is the phase space for four pions of momenta p\, . . . ,p 4 . F 4lT denotes the current-4-7r 
vertex and is the six-pion scattering amplitude. The three-loop diagrams contributing to 
■FjjSi(s) are shown in fig. 1. 

Chiral logarithms are produced by integrations over intermediate momenta with mass- 
dependent boundaries. We find that, in order to calculate the terms proportional to L, we 
can expand the integrand for small M and keep only the relevant terms. It is then crucial to 
find a phase space parametrization which allows us to perform analytical integrations after 
this expansion. We illustrate here a suitable choice for d^ 4 . 

We reduce the 4-body phase space defined as 

dMs; P i, ..,P4) = (2-) 4 * 4 (p - ft iSkp ■ ^ 



8=1 



1 (2vr)3 2 ^ ' 
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(a) 



(b) 



(c) 



(d) 





(e) (f) (g) 

Figure 1. Three- loop cut diagrams contributing to the scalar i^ ne i(s) at order p 6 . 

with s = P 2 to the following product of two-body phase space factors (see for example sec. 
43 of ref. [17]) 



dk 2 dq 2 

d$4(s;pi,p 2 ,P3,P4) = d$ 2 {s;q,P4) X d$2(q 2 ; k,ps) x d$ 2 {k 2 ;pi,p2)-^— — 
k = pi + P2 , q = Pi + P2 + P3 , Pi = p\ = p\ = p\ = M l ■ 



where 

This gives 
d<S>A 



(4.4) 
(4.5) 



1 \WM,M*) \V 2 (q 2 ,k 2 ,M 2 ) W 2 (k 2 ,M 2 ,M 2 ) 



(47T) 



2s 



2q 2 



2k 2 



dq 2 dk 2 



dn k dtt q dn s (4.6) 



in terms of the well-known Kallen function X(x, y, z) = x 2 + y 2 + z 2 — 2{xy + xz + yz). Q s 
is the solid angle spanned by the unit vector q = q/\q\ in the center-of-mass frame of the 
two final pions. Q q is the solid angle spanned by k in the frame where q = and fife is 
the solid angle spanned by p\ in the frame where k = 0. More details on this phase space 
parametrization are given in app. C. The advantage of using this representation is that each 
A function contains M 2 as an argument, which enables us to expand all factors and perform 
all the integrations analytically, at least for the diagrams (a), (b), (c) and (d). 

The integration range for the kinematic variables k 2 and q 2 is determined by the delta 
function which ensures momentum conservation: 



4M 2 < k 2 < 



9M 2 <q 2 <{^~s- M v ) 



(4.7) 



We stress that chiral logarithms stem from both upper and lower integration boundaries due 
to the functional form of the phase space. 

Let us now consider the scalar form factor. For the diagrams (a) to (d) we have been 
able to determine analytically both the values in the chiral limit and the coefficients of L for 
Im F(s). The latter ones are given by 

„2 



Si 



(4ttFY 



(4i 
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where the numerical factors are 

<5a = ^TT, <5 b =— 7T, C =— TT, S d = --TT (4.9) 

o 12 o 4 

For the remaining diagrams the structure is too complicated to perform all integrations ana- 
lytically. In order to extract the coefficient of the leading chiral logarithm from the remaining 
diagrams we set up the following procedure. We performed analytical integrations whenever 
it was possible and the remaining integrations were done numerically. For each diagram we 
then generated a large number of points within a certain range of values of M 2 / s and fitted 
the values of the amplitude minus its chiral limit value with a functional form dictated by 
XPT . Our optimal choice for the number of points, the range in M 2 /s and the truncation 
in powers of M and s of the fit function was determined by our ability to reproduce the 
coefficients Si for the diagrams (a) to (d) within one per cent. 

Summing up the contributions from all seven graphs we obtain the coefficient of L in 
ImFi ne i(s) at three loops: 

s 2 

S - _.. with 5 = (-0.53 ±0.05) tt (4.10) 
[AttF ) 4 

The error on 5 takes into account the uncertainties both in the numerical integration and in 
the fit. The pion mass dependence resulting from the sum of the seven diagrams is definitely 
not compatible with a vanishing coefficient for the chiral logarithm. As a check, we tried 
to fit with a high degree polynomial in M 2 but without the leading chiral log and obtained 
unacceptably high values of x 2 . 

We checked our results also using the four-body phase-space parametrization in ref. [18]. 
This is not useful to do analytical integrations but it provides an important check of our 
numerical routine since it involves different kinematic variables and angles compared to our 
parametrization. 

After performing the dispersive integral in eq. (4.1), our result for ai ne i is 



"inel (S) 



CU?) + - X ( In ^ + iTT 
TT \ S 



S 2 



(4.11) 



(4ttF) 4 

where C(^ 2 ) is a combination of 0(p 8 ) LECs which compensates the //-dependence. Assuming 
that C(/x = 1 GeV) is of natural size, comparing the logarithms from the elastic and inelastic 
part at three loops, we find that the factorization breaking effect is about 5% in the range of 
interest for s. In fig. 2 we plot the relative contribution of the inelastic log to the total three- 
loop log (left panel) and the relative contribution of the inelastic log to the sum of the elastic 
and inelastic log up to three loops (right panel). The LECs are set equal to the values adopted 
in the analysis of ref. [8] . We checked that varying them within their phenomenological ranges 
does not change our conclusions. Comparing the two plots in fig. 2, one can see that for values 
of y/s which are not small compared to A x the three-loop result is not suppressed compared 
to the two- loop and one- loop results, which is a clear signal of breakdown of chiral power 
counting. 
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Figure 2. The relative contribution of the inelastic log. The solid line corresponds to C(/i = 1 GeV) = 
1 in eq. (4.11) while the dashed line to C{ji = 1 GeV) = — 1. The bands corresponds to the uncertainty 
on 8 quoted in eq. (4.10). 

If we go beyond the three-loop level, but still use the chiral counting as a guide, we can 
write cki n ei(s) as the following series 



a small, but non-negligible correction to the factorized one. As long as one remains in the 
chiral regime M 2 C s < A^., factorization of the chiral log does emerge as a property which 
holds to a good approximation. If, on the other hand we abandon the low-energy region, the 
chiral counting is not valid anymore and it becomes impossible to even estimate the relative 
size of the inelastic to the elastic contribution since one should sum the whole series (4.12) 
and the analogous one for the elastic part, before being able to draw any conclusion. In 
fact, if we go to asymptotically large energies, factorization of the chiral log does emerge 
again as an approximate property, but the origin of it changes completely and the coefficient 
of the chiral log gets modified, as we are going to discuss in the next section. Evidently, 
as the energy increases, the inelastic contributions become more and more important and 
factorization does not hold anymore, until one reaches very high energies where it shows up 
again in a very different form. 

We also stress that at 0(p s ) in the form factors there will be contributions of four-pion 
intermediate cuts in the tttt scattering amplitude t^ 8 \ These can be additional sources of 
logarithms. 

5 Chiral logs for asymptotic energies 




n=2 



(4.12) 




We point out that the factorized form of F(s) conjectured in H7r\PT, eq. (1.1), is also not 
consistent with QCD factorization for asymptotically large values of s. In this regime it 



is well-known [9, 19, 20] that the pion electromagnetic form factor can be written as the 
following convolution, 



Ms) 




dxdyT(x,y,s)Mx)My) X [1 + 0(A| CD /s, M%/s)] . 



(5.1) 



Here T is a hard scattering kernel which can be calculated in perturbation theory and cf> n is 
the non-perturbative light-cone distribution amplitude (LCDA) for the pion. In the language 
of Soft-Collinear Effective Theory, this corresponds to the matching of the electromagnetic 
current onto the leading effective operator built with four collinear quark fields, which has a 
non- vanishing overlap with symmetric pion states (i.e. containing only energetic modes) [21]. 

In ref. [10] Chen and Stewart studied the chiral expansion of the vacuum-to-pion matrix 
element of the bilocal quark field operator defining the c/) n (x) in eq. (5.1). They considered 
the tower of local axial twist- 2 operators 0^' a related to moments of the pion LCDA: 



where n is a light-like vector, and proved that for each of these matrix elements the leading 
chiral log is given just by the one in F n . Therefore, according to eq. (5.1) and the result 
by Chen and Stewart, the leading chiral log in Fv(s) does factorize for s 3> Aq CD but the 
coefficient ay in eq. (1.1) is —2, not —1 as predicted by FPrxPT. 

6 Conclusions 

In this paper we have scrutinized the foundations of hard pion chiral perturbation theory. 
We did so on the basis of one example, the pion form factors: in this, as well as in other 
quantities, the leading chiral logs are predicted to factorize in the limit M 2 <C s. For the form 
factors, however, this property has been successfully tested by an explicit two-loop calculation 
in xPT . Our aim was to investigate whether this factorization holds even at higher orders in 
%PT . Since going beyond two loops within a diagrammatic approach in xPT is prohibitive, we 
have based our analysis on dispersion relations and unitarity, and have systematically applied 
the chiral counting to the dispersive representation. This approach lends itself to a recursive 
kind of analysis and allowed us to explain why factorization emerges at the two-loop level. 
Moreover, we have been able to show that for a whole subclass of diagrams this property 
holds to all orders. 

On the other hand, if one considers multipion contributions to the discontinuity of the 
pion form factor (in short: inelastic contributions), one can see that these also generate chiral 
logs. We have calculated the coefficient of the leading chiral log in these diagrams at three 
loops and shown explicitly that these do not factorize. Factorization of the leading chiral logs 
in the hard pion regime is therefore not an exact property. As long as one remains in the 
low-energy regime, but takes the quark masses very small (M 2 <C s) factorization of the chiral 
logs is valid to a good approximation, as our numerical analysis has shown. As the energy 



(ir b \o£ a \0) 



dx (l — 2x) k cp n (x) 



(5.2) 



Jo 
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increases, however, the inelastic contributions appear to gradually become more important, 
until factorization effectively disappears. We conjecture this also because the behaviour of 
the form factors for asymptotically large energies is known. As Chen and Stewart [10] have 
shown, one can analyze the form factor dependence on the light quark masses in the Brodsky- 
Lepage formula: they concluded that the chiral log does indeed factorize for the leading term. 
Nothing is known about subleading terms but there is no reason to assume that the chiral logs 
would factorize for them too. The coefficient of the chiral log in the leading term has been 
given explicitly by Chen and Stewart and it is twice as large as what is found at low energy. At 
intermediate energies some sort of transition between the two regimes must therefore happen, 
and we see no reason why factorization should be valid in this region. 

In the future, we plan to extend this analysis to heavy-light form factors which represent 
one of the most interesting areas of application of H-ztxPT. 
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A 7T7T partial waves beyond tree level 

Consider the dispersive representation of the irir partial waves proposed by Roy [22]: 

2 oo «oo 

4(s) = H(s) + E E / dJKg(8,J)1m1${J) , (A.l) 

where I and I denote isospin and angular momentum respectively and kl(s) is the partial 
wave projection of the subtraction term 

kl(s) = 46! + S -^f(2a° - 5a%) + ±6[S} - • (A.2) 

In order to analyze the possible sources for chiral logarithms, we split eq. (A.l) into the 
contributions from the S- and P-waves, the higher partial waves and the integral from a 
cut-off A 2 to infinity: 



where 



tj(s) = k{(s) + t{(s) SP + 4(s) , (A.3) 

4(s) SP = E E / A ds'Kg'(s,s')lni4(s') , (A.4) 
r=oe'=o j4M - 
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and the so called driving terms d\{s) are given by [23] 



2 oo /»A^ * ^ 



2 oo 



ds'K${s,s')Im4',(s') . (A.5) 



A 2 



The expressions of the kernels can be found in ref. [23]. The subtraction term could in 
principle contain chiral logarithms but the combination of scattering lengths 2a(j — 5a 2 , does 
not have any [2]. As we are going to show, neither the S- and P-wave contribution nor the 
driving terms contain leading chiral logs either. 

In the elastic region unitarity relates the imaginary part of the tttt partial wave to its 
modulus squared: 

Imt I l (s) = a(s)\t I l (s)\ 2 . (A.6) 



The integral over the S- and P- waves may in principle give rise to chiral logarithms: 

,A 2 



SP 



+ 



+ 



4M| 
A 2 

AMI 
A 2 

4M 2 



ds'K$(s,s') 



\ aX s' + ^) V 2 " 8s' M 2 ) + 0(M 4 ) + 0(p 6 



(A.7) 



ds'KH(s,s') 



4 K) 2 (s' 2 ~ 8s'M 2 ) + 0(M 4 ) + 0(p' 
lb 



6n 



ds'K I2 (s,s') 



The partial wave amplitudes needed for the form factors are i[](s) for Fs(s) and for 
Fy(s). Eq. (A.7) involves integrals of the same type of those discussed in sec. 3.1. We find 
that for both partial waves the coefficient of lnM 2 /s is proportional to 



[2a° b° - 5a 2 6 2 ] 



3M 2 



[2(6°) 2 + 27(al) 2 -5(6 2 



2\21 



(Ai 



Although the individual terms in this combination are of order M 2 , they cancel and leave as 
leading contribution something of 0(M 4 ) and therefore beyond the accuracy of this calcula- 
tion. 

According to the expressions of the kernels, terms proportional to L can be generated 
only from Imt(s) up to 0{p A ). Therefore d\(s) cannot produce leading chiral logarithms 
because there Imt(s) starts at 0(p s ). 

Since the tttt scattering amplitude at tree level and zero momentum vanishes linearly in 
the chiral limit, it does not contain terms proportional to L. Therefore 

t {2) {s) = 0{M 2 ) + 0{s) . (A.9) 

From unitarity, by applying the chiral power counting to eq. (A.6), it follows that 

Imi( 4 ) =a(s)\tW\ 2 . (A.10) 
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Hence also the imaginary part of the partial wave to order p 4 contains no term proportional 
to L. Using Roy equations (A. 3), we can get the corresponding partial wave from the 
imaginary part. Since our previous explicit calculation shows that neither the integration nor 
the subtraction term produce any unwanted chiral logarithms, the partial wave to next-to- 
leading order has no terms proportional to sL, hence /5 = in eq. (3.16). By induction one 
may reach the same conclusion to all chiral orders provided one does not consider inelastic 
contributions. This shows that for the partial waves expanded in M 2 <C s: 

n-l 

£( 2 ")(s) = i^ 2n \s) + L ^ ^M 2 ^- 1 -^ + 0{M 2 ) (A.ll) 

A:=0 

the coefficient (3 n -i vanishes to all orders. 



B Inelastic contributions: explicit expressions of the integrands 



We list here the single contributions to the imaginary part of the scalar form factor from the 
diagrams (a) to (g) in fig. (1), namely 



1 



(B.l) 



To express the integrands given here below in terms of the variables of our four-pion phase 
space parametrization in eq. (4.6), we refer the reader to app. C. The momenta of the final 
pions are denoted by p e and pj, see fig. 3, and we define Sij := p% • Pj (notice that s e f = 
s/2-M 2 ): 



Yi 



(a) 



Yi 



5 (25M 2 + 8s) 

24F 6 
1 



(B.2) 



{5 [l05M 4 + 2M 2 (55(si 2 + s i3 + s 23 ) + 20(si 4 + s 24 + s 34 ) - 6s) 



6F e (M 2 -q 2 ) 

+8^Si 2 (6si3 + 3S14 + 6S23 + 3S24 + 4S34 - s) + S2 3 (4si4 + 4S23 + 3S24 + 3S34 - s) 
+S1 3 (3S14 + 6S23 + 4S24 + 3S34-S)+4S 2 2 + 4S 2 3 ) } (B.3) 



Y< 



(c) 



Y, 



6F 6 (M 2 - q 2 ] 



3 M' 

2 

1 



{5(3M 2 + 2(si2 + s 13 + S2 3 ; 



10 M + 5(S12 + S13 + S23) + 3(S14 + S24 + S34 



)-.]} 



(B.4) 



{ [l5M 4 + 20M 2 (si2 + S13 + s 23 ) 



W UF G (M 2 -q 2 ) 2 

+4 (3(s 2 2 + s 2 3 + s 23 ) + 2(si 3 si2 + s 23 si2 + si 3 s 23 ; 

X(14M 2 + 10(S12 + Si 3 + S 23 ) + 6(su + S 24 + S34) + s)} 



(B.5) 



Y, 



AF 6 [M 2 - (pj +p 2 -Pe 



-|2M 2 (5S12 + 3(s e / - Sl e - S 2e - S\f - S 2 f + S34) 
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Y 



+4(si3 + «14 + «23 + S 2 4 ~ S 3e 
+2(si 3 + «14 + «23 + S 2 A ~ S 3e 
+S24 + S 34 + S e f - S 3e - S 4e 
1 



S 4e ~ S 3 f - S 4 /)) + 4 



•S12 (>34 + S e / - Sif - S 2 f 



SAe ~ S 3 f - S4/)) - (Sl e + S 2e )(si3 + Sl4 + «23 

sif - s 2 f - s 3f - Sif) + 11 M 4 | 



(B.6) 



l 15 M 4 + 10M 2 (>12 + 2si 3 + S14 + S23 + S34) 



(/) 2F 6 [Af 2 - (p! + p 3 + Pi ) 2 ] {M 2 - q 2 

+4( Si 3 (si2 + 3si 3 + S 23 ) + SU (S12 + S13 + 3s 23 ) + S34(3S12 + S13 + S23) 



3,)} 



(9) 



x(3M 2 + 2si2 + 2s 2 3 + 2s 24 
1 

~2F 6 (M 2 -q 2 )[M 2 - {pi+p 2 -p e ) 2 
-3stf + 7s 23 + 4s 2 4 - 3s 2e - 3s 2 / + 3s 34 - 4s 3e 



(B.7) 



I - 55 M 6 - 10M 4 (10si2 + 7si 3 + 4su - 3si e 



4s 



3/ 



4s 



4e 



4s 4 / + 3s 



ef) 



-20 M 2 



3S?2 + «12(5S13 + 4s M - si e - 2si/ + 5s 2 3 + 4s 24 - s 2e - 2s 2 / + 2s 34 - 4s 3e 



-4s 3 / - 4s 4e - 4s 4 / + 2s e f) + sf 3 + S23 - s u s le + si e s lf + s u s 2 3 - 2s ie s 2 3 

-Sif S 2 3 - Si e S 2 A + S 2 3 S24 ~ «14 «2e + «1/ S2e ~ 2s 2 3 «2e ~ «24 S2e + Sle »2/ ~ S 23 S 2 f 
+ S 2e S 2 f ~ Sl e S34 + S 23 S 34 - S 2e S 34 + Sle S 3e - S 23 S 3e + S 2e S 3e + S Xe S 3 f - S 23 S 3 f 
+S 2e S 3 f + Sle S 4e - S 23 S 4e + S 2e S 4e + Si e 84, f - S 23 S 4 / + S 2e «4/ ~ Sle S e / + S 2 3 Sef 
-S2e S e / + S13 (si 4 - 2si e - Sif + 2S23 + S24 ~ 2S2e ~ S 2 f + S3 4 - S3 e - S 3 f - S 4e - S 4 / 

+s e f) - 8 3si 2 (2si 3 + 2s i4 - sif + 2s 2 3 + 2s 24 - s 2 / + s 34 - 2s 3e - 2s 3/ - 2s 4e 

-2s 4 / + S e f) + S12(S13 - Sle + S23 ~ S2e)(2si 3 + 2si 4 - Sif + 2S23 + 2S2 4 - S 2 f + S3 4 
-2s 3e - 2s 3 / - 2s 4e - 2s 4 / + S e f) - S13 (2sie + S 2e ) + S13 (- s u (si e + 2s 2e ) 

+s 2e (2s X f - 3s 23 - 2s 24 + 2s 2 / - 2s 34 + s 3e + 2s 3 f + 2s 4e + s 4 / - 2s e /) 
+sie (2si/ - 3s 2 3 - s 24 + 2s 2 / - 2s 34 + 2s 3e + s 3 f + s 4e + 2s 4/ - 2s e f) 

+S23 (sie(2sif - S23 ~ 2S2 4 + 2s 2 / - 2S3 4 + S3 e + 2s 3 / + 2s 4e + S 4 / - 2s e f) 
+s 2e {2sif - 2s 23 - s 24 + 2s 2 / - 2s 34 + 2s 3e + s 3 / + s 4e + 2s 4 / - 2s e /) 

(B. 



-Sl 4 (2sie + S2e 

C Four-pion phase space and related Lorentz transformations 



In order to perform the integral over the four-pion phase space in eq. (4.6) we need to express 
all the dot products p% ■ pj entering the Feynman rules for the cut diagrams in terms of our 
integration variables. We define angles in the following frames: S s , which is the center of 
mass frame (CMF) of the two final pions, S g , which is the CMF of the particles p\, p 2 and p 3 , 
and Sfc which is the CMF of the particles piand p 2 . The labels of the particle momenta are 
explained in fig. 3. The angles and the corresponding reference frames are shown in fig. 4. 
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Figure 3. The kinematic invariants corresponding to our choice of parametrization: s = (pi + p 2 
P3+Pi) 2 , Q 2 = (Pi +P2+P3) 2 and fc 2 = (pi + p 2 ) 2 - 




P3 






k \ 




Figure 4. Our choice of reference frames and angles. 



The latter ones are connected by the following Lorentz transformations 
These are given by 

10 



A, 



1/ 







(C.l) 



(C.2) 
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and by 



A? 



k 2 +q 2 -M 2 
2^/Wyfq 2 




X 1 / 2 (k 2 ,q 2 ,M 2 ) 

2\fW- s fq 2 
X 1 ' 2 (k 2 ,q 2 ,M 2 ) 



X 1 / 2 (k 2 ,q 2 ,M 2 ) . 
* '- sine 



COS (f. 

1,2 , „2_ A/ r2 



X l / 2 (k 2 ,q 2 ,M 2 ) 
2\fWyfq 2 

— sin</> 



COS0 








+q 2 -M 2 



2\/k 2 Jq- 



2s/WJq 2 



- — — Sill Uu =z — 7=^- oin t/y oin 1/ u =z — ^ wo w am t/t. uuout 

\ 2-JWyfq 2 ft 2^/k 2 ^ Y S 2a/PV? / 

with the usual Kallen function A(x, y, z) = + y 2 + z 2 — 2xy — 2xz — 2yz. Now we give 
the momentum vectors in the frames where they have the simplest form. We define the 
four-vectors p\ and p2 in the frame as follows 

/ 1 \ ( 1 \ 



sin sin 



cos <p cos t/fc — sin f/fc 



cos sin f/jt cos C/fe 



, (C3) 



1 


cr(fe 2 ) sin 9 a 
\a{k 2 ) cos 9 a J 



p'i 



1 



<t(&: 2 ) sin a 
\-ct(/c 2 ) cos a / 



with <j(x) = y/l — 4M 2 /x. We give the four- vector p% in the frame E g 

^ g 2_fc 2 +M 2 \ 



P3 



g2_fc 2 +M 2 



X 1 / 2 (q 2 ,k 2 ,M 2 ) 



cos 6k 



lyjq 2 

X 1 / 2 (q 2 ,k 2 ,M 2 ) . n 

— 7= sin t/fc 



and the four- vectors p4, p e and in S s 



/ s-q 2 +M 2 \ 
2y/~s 



Pi 



\ 





A 1 /2( S)9 2 M 2) 





<r(s) sin 9 e cos 4> e 
a(s) sin 9 e sin e 
^ a(s)cos9 e J 



Pf 



\ 



-a(s) cos e sin e/ e 
-cr(s) sin </> e sin # e 
^ —cr(s) cos 6* e y 



(C.4) 



(C.5) 



(C.6) 



After transforming all these momentum vectors onto the same frame, all Lorentz invari- 
ants are written as functions of the seven phase space variables k 2 , q 2 , 9^, 9 a , 9 e , (ft and 

4>e- 
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